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Abstract 



We consider a phenomenological model for the dynamics of Wilson loops in pure SU (N) 
QCD where the expectation value of the loop is the average over an interacting diffusion process 
on the group manifold SU(N). The interaction is provided by an arbitrary potential that 
generates the transition from the Casimir scaling regime into the screening phase of the four- 
dimensional gauge theory. The potential is required to respect the underlying center symmetry 
of the gauge theory, and this predicts screening of arbitrary SU(N) representations to the 
corresponding antisymmetric representations of the same iV-ality. The stable strings before the 
onset of screening are therefore the fc-strings. In the process we find a non-trivial but solvable 
modification of the QCD2 matrix model that involves an arbitrary potential. 
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1 Introduction 

Our current understanding of the quark confinement mechanism is based on the assumption of 
a linearly confining potential with color electric flux between fundamental charges localized in a 
tube connecting them. Wilson loops are the typical gauge invariant operators, i.e. observables, 
one considers: they represent the order parameters which allow to detect the presence of a 
confining phase. It is therefore of extreme importance to study their dynamics. Loop equations 
were already formulated in an elegant way long ago [1] but it turns out to be extremely difficult 
to solve them exactly, since Wilson loops are quite complicated objects themselves being non- 
local and composite (See [2] for a very interesting review of the topic and connections with string 
theories) . 
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Numerical evaluations of Wilson loops are therefore very useful to shed light on the problem. 
Indeed, extensive computations in lattice gauge theory support quark confinement and are able 
to describe quite in detail the behavior of the static quark potential. Recent studies [3, 4, 5, 6, 7] 
have shown, in particular, the presence of an intermediate region where the string tension of 
the linearly confining potential is proportional to the quadratic Casimir operator (of the color 
group); hence the name "Casimir scaling" . However, as soon as the distance between the charges 
increases, at some point screening by gluons becomes energetically favored and Casimir scaling 
breaks down [6, 8]. We just want to recall that screening in general just shows up as a slight 
deviation from Casimir scaling and that one can also consider charges in higher representations 
of the color group rather than the fundamental. For the SU(2) and SU(3) cases, any such 
charge can be screened either to the fundamental or to the trivial representation; for the case 
of SU(N) with N > 3, however, there can be -stable- strings with string tensions different from 
the fundamental one, which are normally referred as /c-strings. 

It would be desirable, of course, to give an exhaustive and more clear explanation of these 
results from lattice gauge theory simulations. It is well-known, in particular, that Casimir scaling 
is exact in two dimensional Yang Mills theories for all kind of loop sizes. In this specific case, 
however, confinement is a perturbative result and de facto the distribution of flux is random, since 
the plaquettes decouple. In addition there is no transition to the screening regime. Some years 
ago, however, Ambjorn et al [9, 10] made, the suggestions that four dimensional QCD, at large 
distances, reduces to two dimensional Yang Mills because of a sort of "dimensional reduction" . 
The proposal was motivated by the physics of spin systems in random magnetic fields which 
displays such a behavior. In QCD then one assumes that a necessary and sufficient condition for 
confinement is a vacuum with random fluxes. The latter play the role of the magnetic field to 
give effectively dimensional reduction, which could then explain Casimir scaling. This proposal, 
as we will comment below, was already tested numerically in [11, 12] and the results of the 
numerical simulations show good agreement. Still, however, one has to explain the screening 
of the charges. Ensembles of Wilson loops, in particular, were considered and their probability 
distribution was computed. 

A phenomenological model along these lines was proposed in [13], where it is assumed that the 
probability distribution precisely undergoes diffusion on the group manifold of the corresponding 
color group. The model is in principle able to explain the transition from Casimir to screening 
regimes and the analysis was carried out for the SU(2) case. In our paper we extend the analysis 
to the SU (N) case and by demanding explicitly center symmetry of our dynamical system we 
can see the emergence of fc-strings as well. We will translate this proposal in terms of a matrix 
model which is a non-trivial perturbation of the 2d Yang-Mills matrix model where area is no 
longer preserved under diffeomorphisms when going from Casimir to the screening regime. At 
the same time we give a more robust mathematical description of the model. 

The organization of the paper is as follows: in Section 2 we review how to relate the Wilson 
loop distributions to diffusion processes. In particular we recall the assumption of dimensional 
reduction. In Section 3 we proceed along the lines of [13] revisiting the free diffusion on the group 
manifold and generalizing it to the SU(N) case. We show the emergence of Casimir scaling. 
In Section 4 we add a drift term discussing briefly some of its properties and map the problem 
to an equivalent quantum mechanical path integral in the presence of a potential. In Section 
5 we demand invariance under center symmetry and see how the model is able to explain the 
transition from Casimir scaling to the screening regime. We also show the emergence of /c-strings 
and some specific choices for the potential. In Section 6 we discuss a candidate effective action 
to describe the Wilson loop dynamics. We consider then the diffusion process at small times 
and make connections with weakly coupled Yang Mills along the lines of [13]. In Section 7 we 
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comment on the vortex density energy which can be predicted from the model and make some 
links with lattice gauge theory (LGT). We end with some conclusions and open problems. In 
the Appendices we collect some group theory definitions and properties we used throughout the 
paper. 

2 Wilson loop distributions as diffusion processes 

We are first going to review some of the ideas underlying the approach we are following. As 
said, indeed, according to [13] one associates the Wilson loop distribution to a diffusion process 
on a group manifold. Similar proposals, however, already appear in the literature, as we are 
now going to revisit briefly in this section. 

Indeed it has already been noticed in [9, 10, 11, 12] that the value to of a Wilson loop 
undergoes Brownian motion as the size of the loop is progressively increased. In addition, each 
such path as a function of the area represents a single random walk and an ensemble of such 
trajectories gives the spectral density p of the values w of the Wilson loop 



Qualitatively, as soon as the size of the Wilson loop becomes larger, p becomes uniform: the 
gauge field is weakly correlated and the distribution of eigenvalues represents disordered config- 
urations. Most notably, confinement is then an outcome of such uniform distribution. 

To carry out such analysis the hypothesis of "dimensional reduction" was explicitly assumed 
and tested. The latter was originally suggested in [9], [10]. More precisely, dimensional reduction 
is a consequence of the assumption of "stochastic confinement", which states that the QCD 
vacuum is made of a disordered color magnetic fluxes. One can then infer that the IR dynamics 
of four dimensional QCD is given by two dimensional QCD with quite good approximation. For 
another derivation of dimensional reduction see [14]. Dimensional reduction is a consequence, in 
any case, of the specific stochastic assumption on the nature of the QCD vacuum. Note also that 
stochastic confinement implies Casimir scaling, which was conjectured long ago before recent 
lattice simulations. 

Dimensional reduction was then tested in [9, 10, 11, 12]. Qualitatively, at weak coupling (up 
to usual roughening transition problems), dimensional reduction seems to be reproduced. At 
strong coupling, however, a more quantitative and remarkable result was obtained 



namely the spectral density p of a Wilson loop of size / x J in four dimensions at coupling (3 
depending on the value of the eigenvalue a (consider here for simplicity the SU (2) case where 
there is only one eigenvalue as discussed in more in detail below) is roughly equal to the two 
dimensional one evaluated at a coupling (3 2 d which is a function of the four dimensional one 
f3. Note that the two spectral densities almost coincide already with the renormalization of 
one single parameter. Summarizing: one is studying a two dimensional effective model where 
the "residua" of the four dimensional asymptotic freedom are in the dependence of the two 
dimensional coupling (3 2 d from the four dimensional one /3. In the approach of [13] this effective 
model is described by a diffusion process on the group manifold and the dependence from the 
coupling is reabsorbed in the "time" t entering the diffusion equation. 

The proposal of [13], in particular, is to identify the spectral density p with the kernel of 
the the heat equation on the group manifold of the color group. Computations have been 



p(w) = (0 | S(W -w)\0) 



(1) 



P/x5(«>/3) ~P/xj(«>Am) 



(2) 
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performed assuming somehow the dimensional reduction already at work, since they turn out 
to be computable in a deformed version of QCD2, where character expansions on the group 
manifold are exact. Actually it is also well known that the partition function of 2d Yang-Mills 
satisfies a heat equation on the group manifold so the proposal to identify the spectral density 
p with the heat kernel solution on the group manifold of the color group is not surprising, once 
dimensional reduction has been invoked. The same conclusion is obtained from the geometrical 
lattice gauge theory (LGT) actions discussed in Section 7 where we believe however it comes 
out in a more transparent way. We want to stress, however, that our theory is not QCD2 but a 
perturbation of it, and the non-trivial physics lies in this perturbation. 

3 Free diffusion of Wilson loops on SU (N) and Casimir scaling 

We are first going to review the free diffusion of the spectral density of Wilson loops p for the 
SU(2) case which was already discussed in [13]. We then generalize to SU(N). Free diffusion 
will give Casimir a scaling regime, i.e. the string tension of the confining linear potential is 
proportional to the quadratic Casimir of the color group. 

Casimir scaling in general holds, approximately, at intermediate ranges, from onset of con- 
finement till color screening 1 . Theoretical arguments in favor of it are dimensional reduction 
as discussed and also Witten analysis in the large N limit [15], where Casimir scaling becomes 
exact. Numerical evidence was provided quite recently by lattice simulations [4, 5]. 

3.1 Revisiting free diffusion on SU(2) 

This is the case discussed in [13]. It is useful to review it and make some additional remarks 
before moving to the more general SU(N) case. 

The probability distribution of Wilson loops is given in this case by 2 

P(9, t) = sin 2 6 G(9, t) (3) 

where G(6,t) is the kernel of the heat equation on SU(2) = S 3 , i.e. it satisfies 

(|- A ) 0( "' () = ^ W(t> - (4) 

In this case 9 is the first polar angle on S 3 while t is time. The Laplacian on S 3 is: 

1 sin 2 9 1 

A s 3 = n ■ 2/i d e ( sin2 e d e) + o ■ 1 d 4> ( sin 4>d ( j > ) + - sin 2 9 sin 2 <f> d 2 , (5) 
2snr0 7 2sm<p 2 v 

where the factor of ^ is conventional. There is an SO (4) symmetry group acting on the S 3 , 
and the eigenfunctions of this Laplacian are the y n im(9, </>, </?), / = 1, . . . , n, m = —I, . . . , I, with 
eigenvalues 

A 5 3^nZm = - 2 n ( n + 2 ) ynlm , (6) 

1 Roughly till 1.2fm [5]. 

2 In our case, G is the spectral density p. Compare with [12]: their p = -|G. In addition our time t will depend 
both on the coupling and the size of the loop as a consequence of the model we are considering. We give some 
explicit examples of this dependence in Section 6. 
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and n = 0, 1, 2, . . . The Cartan subalgebra of SU(2) however only depends on a maximal torus, 
parametrized by 9. Assuming gauge invariance of the background and of the Wilson loop 
operators fixes the remaining angles 4> and (p to constant values. Thus, we will be considering 
^-dependent functions only. Instead of using the it will then be more convenient to use 

a basis SU{2) characters, which are eigenfunctions of the "radial" part of the Laplacian: 

A s s Xj (e) = -2j(j + i) Xj (e) (7) 

where n = 2j. The SU(2) character is: 

The same will go through for SU(N). Again one can invoke gauge invariance arguments. 
This restricts then the dependence of G only on N — 1 "radial" coordinates of the Cartan 
subalgebra. We know that from the point of view of LGT, as we will see in section 7, this is 
natural, since any function one builds up out of group elements has to be a class function in order 
to preserve gauge invariance, i.e. it has to depend only on radial coordinates. This reduction 
to "flat space" namely the maximal torus spanned by the N — 1 Cartan generators has been 
already pointed out in [16] where quantum mechanics on group manifolds was discussed. There 
the reduction to radial coordinates can be interpreted geometrically as rotational invariance in 
group space sense. In more mathematical language this is related to the so called intertwining 
operators [17]. 

Using standard spectral decomposition methods the kernel G is easily found to be 

Glf,t) = \m± „^ e -<»^. (9) 

n=l 

We have shifted n by one unit, so n = 2j + 1 and we see that the exponent indeed falls off with 
the Casimir. We now compute Wilson loops in this Brownian motion ensemble average. One 
easily finds finds 



log^) hih + l) 



log(W j2 ) J2O2 + I) 



(10) 



where (Wj) is the vev of the Wilson loop operator in the spin j representation of SU(2). One 
sees therefore that by definition we are in the Casimir scaling regime. Note also that as t — > 00 
G — > 2/ it, i.e. one has a uniform distribution which is an indication of confinement as recalled 
in Section 2. 



3.2 Free diffusion on SU{N) 

We now generalize the diffusion process to the SU(N) case. As explained above, gauge invariance 
allows us to restrict the Laplacian to the Cartan subalgebra, and we get the following free heat 
equation: 

(i- A ) 0(t ' 9) = 7^) s<(W) - 

Here = 0i = (6\, ...,6n) are the Cartan coordinates, subject to the SU(N) constraint: 

N 

J> = 0; (12) 
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in the following we will always assume this constraint and also drop the subscript from 6*j. J(9) 
is the Jacobian associated with SU(N) Haar measure. The explicit expressions and more details 
of the solution are given in Appendix A. 
Let us consider the more general problem 

(Jt - a) g(t; 9, 6') = 5(t)5(e - 0') . (13) 

The right-hand side gives us the choice of boundary condition at t = 0. The solution therefore 
is: 

g(t;e,9') = 9(t)J2xx(0)xx(0')e^ C2Wt . (14) 

A 

where the sum runs over all irreducible SU(N) representations, labeled by A. In turn, this is a 
sum over integers, as in the simple SU{2) case (9). Details are explained in Appendix B. 

The above is obviously the partition function of 2d Yang-Mills on the cylinder, where t = 
S'ldYM^ where A is the area. Sending 9 = 9' — > 0, we get the partition function of QCD2 on 
the sphere. Both theories are well-known to have phase transitions in the large N limit. In this 
paper we will be concerned with the case 9' = 0. In that case, the above correlator corresponds 
to 2d Yang-Mills on the disc: 

g(t,9) = 9(t)^2xx(0)dlm(X)e-^ c ^ t . (15) 

A 

g(t, 9) can also be seen as a Brownian motion probability p(t, 9): 

g(9,t)=p(9,t) . (16) 

This relation will get modified in the interacting case. 

We will insert a Wilson loop operator with holonomy specified along the boundary of the 
disc, so the picture where we identify the partition function of 2d Yang-Mills with the Casimir 
dynamics dimensionally reduced to the surface spanned by the loop is correct in the free case. 
In the interacting case to be discussed in the next sections this picture will essentially remain 
unchanged; only the theory on this surface will no longer be pure 2d Yang-Mills. 

The Wilson loop average is now given by computing the expectation value of the Wilson 
loop operator in the ensemble (16): 

(W x (9)) = J[dO]W x (0)p(t,9) . (17) 

Here and in what follows d9 will be understood to be the Haar measure on SU(N): [d9] = 
d N ~ 1 9J 2 (9). Let us now discuss what the Wilson loop operator should be. Clearly, gauge 
invariance requires it to be a class function of 9. That means that we can expand it in a 
basis of SU(N) characters. Further, the holonomy around the loop singles out a particular 
representation. Therefore, up to an overall normalization constant we get 

W x (9) = X \{0) • (18) 

This is of course a usual fact in 2d Yang-Mills. However, notice that the above argument applies 
to the interacting theory as well. We can now easily compute its expectation value in the free 
case: 

(W x (9)) = dim(A)e-5 c '2(A)t ? (19) 
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where we used standard properties of characters, summarized in Appendix B. Notice that this 
directly leads to Casimir scaling as in the SU (2) case, (10). The factor of dim(A) is an irrelevant 
normalization constant, which depends on our overall choice of normalization in (18). The 
crucial fact here really is that the probability p(t, 9) is correctly normalized, 

J[d0]p(t,0) = l, (20) 

so that the ensemble (16) is a good Brownian motion ensemble. 

All of the above is completely standard in 2d Yang-Mills. When we include the potential, we 
will see that the representation of the holonomy of the Wilson loop will no longer be the only 
one contributing to its expectation value. There will be a mixing where all representations of the 
same 7V-ality will contribute, and this is of course what makes the dynamics of the interacting 
theory non-trivial. 

4 Diffusion with a potential and screening 
4.1 Drift picture 

As we have seen above, free diffusion on SU(N) reproduces Casimir scaling. In QCD, however, 
the external charges interact with intermediate strings that can be in arbitrary representations. 
This is reflected by the fact that Wilson loop averages are no longer of the type (19). We expect 
states in arbitrary representations to contribute to the final answer for the expectation value. 
Thus, we need to generalize the free diffusion equation (11). An interaction term [13], namely a 
small drift term V(9), will give us a Fokker-Planck equation that is no longer diagonal: 

^p(t, 9) = Ap(t, 9) + d t {d*V P (t, 9)) . (21) 

The fact that the interaction is given by a drift term ensures that there is a conserved current 
and the heat equation takes the form of a continuity equation. This is standard and we will not 
dwell on it. 

We explained above that the dynamics of the theory mixes the representations that contribute 
to a Wilson loop. The contributing representations are not arbitrary though; there is a restriction 
coming from their A^-ality. iV-ality is the number of boxes of the Young tableau modulo N, see 
Appendix B for the precise formula. For example, it is well-known that for large loops the 
leading contribution is that of lowest energy, which is the antisymmetric representation of fixed 
A-ality. Thus the behavior is quite different from (19). We will see that this A^-ality dependence 
follows from a drift that respects the center symmetry of the background. 

Let us briefly review the SU(2) case discussed in [13] before we do the SU(N) generalization. 
The diffusion equation with drift term can be rewritten as: 

(i- A ) p<t - < " = ^ s(tW " <22> 

where we incorporated the boundary condition as a source on the right-hand side, and 

A = A + AV + d l Vdi . (23) 
is the modified Laplacian. Note that this operator is no longer hermitian. 
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To solve this equation we proceed as before; we expand in a basis of eigenfunctions ip n like 
in (14): 

p(t, 9) = 9{t) Me)M0)e- E " t+v( V (24) 

n 

where the constant term was added for convenience. The new solutions ip n are eigenfunctions 
of the modified Laplacian: 

A^ n = -E n ^ n , (25) 

and the boundary condition requires 

J2 M9)M0)e v ^ = - Jj- 8(6) . (26) 



The eigenvalue equation (25) amounts to 



1 d { d 

sin 2 9—J) n + sin 2 6—ip n )= -E n ^ n . (27) 



sin 2 9 89 \ 89 Y 89 
It has an obvious zero energy solution: 

^ = e -vW . (28) 

Notice therefore that the asymptotic distribution in the diffusion process is no longer uniform 
on S 3 as in the free case. 

As discussed, we now require that the drift respects center symmetry (we will explain the 
role of center symmetry in more detail in section 5). We show in Appendix B that the center 
acts on 9 by a shift 9 — ► 9 + ir when N = 2. The key fact here is that the Wilson loop operator 
transforms under center symmetry according to its iV-ality: 

W j (e + ir) = (-) 2 'W j (9) . (29) 

Loops in integer and in half-integer representations will thus be screened differently. In general, 
however, the drift term does mix the representations. If we want it not to mix them, it must 
have the same symmetries as the background: 

V(9 + tt) = V(9) . (30) 

This is nothing than reflection symmetry around the equator 3 . 

The above symmetry gives at asymptotic times (i.e. strong coupling) screening of the integer 
Wilson loops to adjoint loops and screening of half integer loops to fundamental loops [13]. The 
role of this symmetry will be simply to keep, at asymptotic times, states with different iV-ality 
distinct. Note however that the presence of the drift by itself does not ensure a perimeter 
law. In the SU (2) case, one expects a perimeter law for the integer representations (this would 
correspond indeed to a k = string). Still one has area scaling as can be checked. To implement 

3 This of course reflects the fact that SU(2) is the double cover of SO(3). In the usual SO(S) representation 
where ip is the rotation angle around the azymuthal axis of the two-sphere, 6 — i (p. ip has periodicity 2tt, which 
is why 6 only has periodicity 7r. 
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the perimeter scaling one should generalize and require an explicit time dependence of the drift. 
In that case, asymptotically, one could impose [13]: 

V(e,t)^V 00 (9)e- m ^ t (31) 

with m mass of the "gluelumps" . 

Before moving to the SU(N) case it will be useful to map the diffusion problem in the 
presence of drift to a diffusion process in the presence of a potential, to be thought of in path 
integral language. It is well known that a Markov process with some specified drift can be 
mapped to a quantum system in path integral theory. Here we simply follow the same route 
working however in the Euclidean at the level of the quantum system. This mapping will be 
explicitly used in the next section. We first review and give a slight generalization of the ST/ (2) 
case discussed in [13]. 



4.2 The screening potential 

The transformation that maps the diffusion problem with a drift to the quantum mechanical 
problem of a particle in a potential is easy to find: 

MO) = e-^M9) . (32) 
It is easy to check that the diffusion equation 

(| - A) p(t; 6, 6>) = -^6(0- 9>)6(t) (33) 

becomes 

(|"A + U(6)j g(t', t; 6', 0) = -^~ e 5(6 - 9')6(t - t') 
where the potential is given in terms of the drift as: 



(34) 



U{e) = ~AV + - A {V'f = e -H») Ae^W . (35) 
The rescaling of the wave-functions (32) also induced a rescaling of p(t, 9) according to: 

p(t,9) = g(0,t;0,9)e-^ v ^ v ^ . (36) 

Using (32), we get: 

g(t', t; 9', 9) = 9{t - t') £ M9)M&) ■ (37) 

n 

This also explains our choice of overall scaling constant in (45). 

The above analysis straightforwardly generalizes to the SU(N) case. The heat equation with 
drift (22) can again be solved in terms of new eigenfunctions 

p(t, 9) = 9(t) M0)M0)e- E > t+V ^ . (38) 

A 

The sum ranges over all irreducible SU (N) representations, which is again a sum over integers 
as explained in Appendix B. 
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(40) 



As before, there is a non-uniform zero mode 

^ = e - y W . (39) 

By the same mapping (36), we get: 

(jt - A + 17(0)) </(t, i'; 0, 0) = 5(9 - <f>)8(t ~ , 
with the potential given by 

U(9) = e-^^Ae^W . (41) 

Again, (40) is solved by 

g(t,t'-0,ct>) = O(t-t')Y,M())M<t>)e Exit - t ' } , (42) 

A 

where tp\ and ip\ are related by the same wave- function rescaling (32). The zero mode is 

■fo{0) = e~v v W . (43) 
Plugging (42) back into (40), we get the conditions on the eigenfunctions: 

(A-U(6))M0) = -ExM0) • (44) 

Our task in what follows will be to solve this equation, subject to the boundary condition that 
follows from (40): 

J2M0)M0') = ^S(9-9'). (45) 

A J 

In the next section we show that the quantum mechanical model (44) is exactly solvable, and 
explain how the solution is connected to the physics of string iV-ality. 

5 Center symmetry the emergence of fc-strings 

We have seen that the presence of a small drift term and/or equivalently of a potential can bring 
about small deviations from Casimir scaling which are expected when screening of the sources 
takes place. We are going to develop this point further and show explicitly the emergence of the 
fc-strings. The only assumption one has to make is invariance of the potential and/or drift term 
under center symmetry. 

To that end, let us briefly recall the role of the 7V-ality of the center of the color group. In 
our case the color group is SU(N) and the center is Zjy. An element of will act on a SU(N) 
representation via a factor 

exp {—^j ■ (46) 

We call k the 7V-ality of the representation. Now it turns out that fields which can give string 
breaking are the ones with iV-ality different from zero. Therefore confinement can be stated as 
follows: start from pure gauge SU(N) theory plus matter and consider those fields with iV-ality 
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different from zero. Send their masses to infinity. Confinement amounts to show that in this 
limit there exists a confining phase in which the work needed to separate a quark and an anti- 
quark to a distance L tends to oL where now the string tension a depends on the representation. 
In the screening phase, in particular, that is to say after string breaking, it will depend only on 
the N ality of the representation. 

In general, matter with non zero iV-ality breaks the global center symmetry of the lagrangian 
of the theory, which is required by the standard order parameters (Wilson loops, Polyakov loops, 
't Hooft loops and similar) to be able to detect the confined phase. This is why the assumption 
of invariance under center symmetry is natural as well in the case of the diffusion process which 
is supposed to describe different regimes in the confined phase. 

5.1 String tension from the diffusion process 

We now come to the core of the problem, which is the SU(N) generalization of the interacting 
diffusion process. 

We should first discuss how the center of the gauge group acts on the Wilson loop operator. 
In the SU(N) case, this operator is the SU(N) character in representation A, (18). The center 
acts on it as: 

W x (0i + ^) = e™ h W N W x {9i) , (47) 

where k(X) is the iV-ality of the representation. The proof of this is given in the Appendix B. 
The background and the potential, however, are invariant under center symmetry, as already 
discussed. This means that the drift satisfies: 

v{e i ) = v[e i + 2 ^ . (48) 

Hence, since the background and the potential are invariant under center symmetry, the above 
property (18) is a fundamental property of the loop that distinguishes states of different iV-ality. 
More precisely, we find asymptotically: 

(W x (0i))^Ce- E ^ , (49) 

with some proportionality constant C which will depend on the representations. E„ is the 
tension of the string. For example, for SU(2), = for A = j G Z, whereas it is non-zero if 
j € tjZ. In other words, loops in half-integer representations are screened to fundamental loops, 
and loops in integer representations are screened to adjoint loops. In the SU(N) case we expect 
the same behavior (49), where now representations can differ by their iV-ality. Concretely, what 
we expect is that is the string tension of a string in representation ft, which is the totally 
antisymmetric representation of the same N- ality as A. 

5.2 iV-ality of the eigenfuctions of the diffusion process 

To find the expectation value of the Wilson loop, we proceed as follows. As before, the expec- 
tation value if computed by integrating the loop operator over the diffusion ensemble 

(W x (9)) = J[d9]W x (9)p(t,0) , (50) 
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where we are using the Haar integration measure (see Appendix B). p(9; t) is now the modified 
diffusion probability (38) expressed in the new set of eigenfunctions ip\ satisfying the Laplace 
equation (44) with the potential: 

(A-U(9))M0) = -ExM0) ■ (51) 

In this picture, the E\s are obviously positive definite because U{9) is (see Appendix C). From 
now on we will often supress the ^-dependence, but this should be understood in all the character 
formulas. We also need to express the Wilson loop itself in the new basis. Thus, we define 

W x (9) = X x(0) 

= 5>A^W- (52) 

Finding the is now equivalent to finding the coefficients u^. 

There is a restriction on the eigenfunctions in which we expand the Wilson loop in (52). 
Remember that the Wilson loop transforms as (47) under center action. It is not hard to see 
that the eigenfunctions must respect this property as well: 

M0i + jj) = e-»~M0i) ■ (53) 

This also means that the representations we sum over in (52) all have the same 7V-ality, in other 
words we sum over Young diagrams with a fixed number of boxes modulo N. Thus we impose 

w£ = if k^K . (54) 

Since this property is built into the definition of we will not indicate this restriction explicitly 
in the summation signs below. 



5.3 Solution of the diffusion equation: from Casimir to screening regime 

We now explicitly solve the Laplace equation (51). First of all we invert (52): 

n 

Obviously, c is nothing but the matrix inverse of u: 

£« = 5>^ = <5£. (56) 

V V 

Like before, vanishes whenever k^^ k v . 

It is now useful to expand the wave equation (51) into modes. Explicitly: 

= ^c^N^Xu (57) 

fiua 
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where we used 

X,(0)xx(0) = Y J %xXu(9) . (58) 

V 

The iV's are the tensor multiplicity coefficients for SU(N), also called the Littlewood- Richardson 
(for details, see Appendix C). Remember that U(9) has to preserve iV-ality; therefore, 

u c = if k a ^0 . (59) 

Filling the above in (51), we get the following general equation: 

fl (TV 

Notice that this has now become a purely algebraic equation. We can get a closed form for 
the coefficients by multiplying the whole equation with (-,Xn) an d using orthogonality of the 
characters. We get: 

c^A-CaMHEE^^S* ( 61 ) 

The Littlewood-Richardson coefficients preserve iV-ality, that is, N£ v = if / k a + k v . 
Therefore, since k a = and c also preserves iV-ality, the above equation is trivial unless k^ = k\. 
Taking fj, = A in (61) and multiplying with u, we get an explicit expression for the energy: 

E x = ^y x Uc 2 (p) + ^u a NA . (62) 

liv \ a J 

At this point it is very convenient to introduce the following notation: 

i£ = 2>*W£, (63) 



cr 



(for its properties, see Appendix C). Since N% v preserves iV-ality, it is clear that n preserves it 
as well: 

< = if k^k v . (64) 
We define also the following bilinear form: 

< = ^C 2 (^)+<, (65) 
which by definition is iV-ality preserving. The energy can then simply be written as: 

Going back to the full equation (61), we can rewrite it as 

E X % = (c-'gcK (67) 
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and filling in the solution we found for the energy, 

{gcY x = (c- V)M • (68) 

This is an eigenvalue equation. Notice that the index A labels the eigenvalues. Thus, regarding 
c as a vector, = c^, we get the eigenvalue problem 

9 -C(A) = (C~V)(A)C( A ) (69) 



where we now regard the c's as eigenvectors of g. Thus, we have reduced the problem of solving a 
differential equation on the group, to the algebraic problem of finding eigenvalues of eigenvectors 
of g, which is a given matrix once we choose the potential, see (65). For simple potentials, g is 
almost diagonal and this equation can be readily solved. 

We saw that the ground state of the system has zero energy. This implies: 

m = m Mm ■ <70> 

Combining this with the expansion for U (57), after some manipulations we get: 

CqU P = ^2 c Qu u u = , (71) 

V V 

where we used the charge conjugation properties of c derived in Appendix C. 
5.4 Linear fluctuations 

For the purposes of this paper, it is enough to solve the diffusion equation explicitly at the 
linearized level. Thus we can linearize equation (61). We set: 

< = K + z< 

Ex = C 2 (X) + eSE x 

u<r -»■ u' a = tu a (72) 

and work to first order in e. Obviously, to lowest order in e the equation is trivially satisfied: 
this is the unperturbed solution. At the next order, we get two equations: 

6E X = J2 U ° N ^ 

a 

5E x 5» + 5c»(C 2 (\)-C 2 (v)) = E^a- (73) 



This is easily solved: 



^ = 5W^M?"'^ (1 -«' (74) 



k„=o 



In other words, if hy is the non-diagonal part of riu, 



c ^ = AM I 

* = V + cJicM *' (75) 
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Notice that the 5cp are not determined from the above. This is reflected by the fact that the 
fluctuation is expressed in hy not riy. This is due to the fact that we can always change them 
by appropriate normalization. In particular, we can choose b~c\ = 0, in other words c\ = 1 at 
least to second order in perturbation theory. 

Boundary conditions 

We still need to discuss the choice of boundary condition (45). It can be rewritten as: 

£ Xx(9) dim(A) = £ Va Wa(0) . (76) 

A A 

Working this out, we get 

dim(/,) = ^<dimH . (77) 

Taking ji = 0, this reduces to 

£c5tadim(i/) = l. (78) 

The full boundary condition equation can be rewritten as 

£( C J-c£)dim(A) = 0. (79) 

A 

We can now analyze these boundary conditions at the linear level. We get: 

£ CM ~ c 2 (x) ( dim(A) < - dim( ^ <) = • (80) 

A 

When |U = 0, we get the following condition on the potential: 

^C^A) (dim(A) " 1)uA = ' (81) 

A 



5.5 Extracting the fc-strings 

We have now solved the model in the previous section in implicit form, the solution being given 
in equations (66) and (69). For linear perturbations, we obtained a closed form for the solution 
in (73) and (75). We are now ready to compute the expectation value of the Wilson loop and 
extract the tension of the minimal string. 

Before going into details, let us briefly describe what will happen. The important point 
is the following: at large t, because of the exponential damping, the leading contribution will 
be the one with lowest energy E^. We are summing over all representations /j, that have the 
same iV-ality k\ as the loop operator W\. The leading contribution is therefore given by the 
representation /j, that minimizes E^. Since the leading term is 

E lt = C 2 (i*)+0(e) , (82) 

we see that the leading contribution to a Wilson loop (W\) will be given by the representation of 
lowest Casimir. But that is the totally antisymmetric representation with iV-ality k\. Thus we 
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have found the /c-string! Thus, we get that a Wilson loop in representation A is screeened to a 
loop in representation [i, where is the totally antisymmetric representation with k\ quarks. Of 
course in the SU(2) case we reproduce the result of [13], but we see that in general this model is 
very powerful in predicting the screening behavior of the background. This conclusion remains 
unchanged taking into account small perturbations in (82), and in general any perturbation 
produced by the potential, as long it keeps the ordering of the lowest-lying energy levels. 

We now give the details of the above. All the ingredients to compute the Wilson loop average 
have been given earlier: 

(w x (e)) = |[d0]w A (0) P M) = ^M^^ 
Vy(o) 

V'o(O) 



Notice that in this computation it was important to use the fact that p(t,6), rather than the 
correlation function g(t,9), is the correct diffusion density; one indeed easily checks that it is 
normalized to one, 

J [d9] p(t, 9) = |H e~ Ext J m M0)MO) = 1 • (84) 

A 

Let us now look at what happens for large t: 

(^W) = ^^e-^ (85) 

where A(A) means the antisymmetric representation of A-ality k\. In the case that A is in the 
same A-ality class as the trivial representation, the loop gets screened to: 

(W x {9))=u° x , (86) 

that is, a constant. 

The above computations were exact. We now work out the linearized case explicitly. We 

get: 

9(f, 9, 9') = £ Xx(0) Xx (9')e-^W^ x )t + £ Xa Wxa^O ( .y e _c 2 (A). _ h a e -c 2 (x> )t) (g?) 
^ j-^ 2 (A) - G 2 (A ) 

The additional term correctly respects the symmetry in 9 and 9', as it should. 
5.6 Specific choices of the potential: examples 

The approach of [13] is phenomenological, and extra input is required to determine a phenomeno- 
logically intersting potential. For illustration purposes, and since the model with the potential 
is quite interesting in its own right, we now give some examples of simple potentials. A first 
rather trivial example is where u is only non- vanishing for the trivial representation: 



u. 



= u5 n 
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where u is a number. We get 



Ex = C 2 {X) + u. (89) 

Thus in this case we see that c remains unmodified, and the whole effect of the potential is to 
shift the zero-point energy of the theory. From the point of view of 2d Yang-Mills, this is adding 
a cosmological constant term to the action. 

We now look at the following less trivial example: 

U a = UXa{4>) , (90) 

the SU (N) character for some fixed 4>, that is, a collection of numbers, and u is again an overall 
constant. There are two limiting cases: 

Xm(°) = dim(/i) 

Xn(p) = dim 9 (/i) (91) 

where in the second case we have conjugated the group element <\> to the trivial representation 
p by exponentiation with q. We get: 

n<t = u X MxM) • (92) 
This result is exact. Linearizing, we get: 

SE X = u( x> ,(<P)) 2 

< ~ - cm-cm **^- (93) 

where /i/y. The expectation value of the Wilson loop is now: 

(W , ( 0)) = e -E- x t _ u xM. v ^(^)xm(o) C2M t m 

Notice that the correction to the energy in this simple example is the square of the quantum 
dimension of the representation A. Asymptotically this leads to the quantum dimension of the 
antisymmetric representation of iV-ality A, which is somewhat reminiscent of the sine law. We 
will get back to this point in the conclusions. 



6 Wilson loop effective action and small time limit 

6.1 A proposal for an effective Wilson loop action 

As recalled in the introduction it is quite hard to study the Wilson loop dynamics. One of 
the outcomes of the approach we are following is that one can write an effective action for the 
Wilson loop dynamics. The SU{2) case was already treated in [13]. Here we revisit the problem 
and sketch the SU(N) generalization. 

Before doing this we would like to recall that the mapping from the description with a drift 
term to the one with the potential can be obtained in a compact and elegant way as follows: 
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start from the the operator A in (23). To go to a hamiltonian description with a potential 
perform a similarity transformation 

H = TAT' 1 (95) 

with the operator T given by T = exp(V/2). One can easily check that the hamiltonian H will 
be given by A — U, with U precisely given by (35). Note that (up to normalization constants) 
if one defines 

C 6i = d 6i + \de i V; C+ = -d 6i + \d 6i V (96) 
then the hamiltonian can be written as 

i 

This allows to prove more easily that H is a self-adjoint operator with real eigenvalues. In 
addition the ground state xo is given by Cg t xo = which once solved gives precisely exp(— V/2). 
This is the corresponding infinite time (i.e. equilibrium state) in the Fokker-Planck language 
[18]. 

Consider now the effective action. In the SU{2) case the [13] proposal was to write the 
probability distribution as 



r 6(t)=6 

P(lf,t,6',9)= [dO(t)} e - 5 ef#(*)l . (98) 

J0(t')=0' 



We went back to the probability distribution defined in (3). This is related to the two-point 
function g(t',t,6',9) by the Jacobian. Tthe effective action is given by 

S cS [9(t)} =Jdt (^9 2 {t) + U{9{t)) \ . (99) 

Let us review briefly before moving to SU{N) some notions from the theory of stochastic 
differential equations [18]. Start from Langevin equation 4 

e\ = -\h{e{t)) + m(t) (ioo) 

where fi is a function of 9(t) and r\i the noise (which we assume Gaussian). Then the probability 
g(9,t) which can be show [18] to satisfy the Fokker-Planck equation can be expressed as 

N-l 

P(9,t) = (H5[9 t (t)-9 t ]) v , (101) 

i 

where the average is taken in this case with respect to the gaussian noise. 

4 Here we assume i = 1, .., iV — 1 as in our set up and the 9i(t) are then the coordinates in the Langevin 
description. They have -no- relation with the 6i (time independent, i.e. just coordinates) appearing in the Fokker- 
Planck equation we have been considering. The coordinates in the effective action (99) (which are time dependent 
indeed) are thus the one in the Langevin description. We will always put explicitly the time dependence when 
we consider the 9i(t) entering in the Langevin description. 
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The effective action appearing in (98) is given by: 

/•*" 1.1 1 

Sesim] = J dt 2 m(t) + ^h{0it)) 2 - n-difMt))} . (102) 

This action, once one assumes fi = QdiV where (f2 is the diffusion coefficient which appears 
in the two point function of the gaussian noise), i.e. a purely dissipative Langevin equation 
which will produce a drift term of the kind we are using, can be checked to correspond to the 
action (99) proposed in [13] for SU(2) where of course we keep only one 9. 

On a general Riemannian manifold, however, the probability 5 

,0(t")=fl" 

P{9",t";9',t') = / d9(t)T\[det(e)]- 1 exp(-S cS ) (103) 

Jo(t>)=e> V 

where now the action has a complicated expression [18] containing vielbeins and their deriva- 
tives. It would be interesting to generalize the study to SU (N) case. 

In addition in the computation above we always set the diffusion coefficient = 1. But in 
this effective action approach it plays the role of h. Restoring it means that some terms will 
be dominant with respect to others. For instance the second term on the right-hand side of 
(102) can be shown to be subleading. This is relevant if one has in mind to do saddle points 
approximations around minima of the potential. If one kept indeed only the leading order terms 
one would not be able to resolve as usual the degeneracy of the extrema. 



6.2 Small time limit: diffusion on the Lie algebra 

Consider the quantity r(p), namely the variation of the size of the loop p with respect to time 
i. 6 To get some more explicit information about this function, we will follow [13] and look at 
two particular limits of the diffusion probability. From the kinetic term we get 

dpr(p)(^) 2 , (104) 



dp J 



where 



rM = (£) ■ (105) 

So we look in the regime where the potential becomes irrelevant and we can use the free diffusion 
problem: 

p(t, 9) = Xx(9) dim(A) e-^W* . (106) 

A 

Since we are considering Brownian motion on the circle [19, 31], we can use Poisson resummation 
to rewrite this as [21] 

(27rt) — a - h J \Pi + ^ii) 



5 Here Stratanovich convention is assumed for the operator ordering. 

6 The size of the loop p has not to be confused of course with the spectral density p{9, t) which we have been 
discussing and in our case is given by G(0,t). We keep this notation to make contact with [13]. 
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See Appendix A for full details and the definitions of D and J. 
In the SU{2) case, this reads explicitly: 



oo 



v (t 0) = - V 6 + 27rn e -{e+^nf/t (W8) 

v ' n=— oo 

It is now possible to take the limit t — > directly. The only term contributing gives: 

which is the result in [13], and we are also taking the 0's to be small. In doing so, we have 
broken the periodicity of 0. This method immediately generalizes to SU{N). Taking t small in 
(107), we get 



P(t,0) * — L_ ^e"^^ = —L^e-(i-i) ^ . (110) 
(2vrt)— (2vrt) — 



D(0i) _X-N fl2/9 , 1 _( 1 N_\ fft. 

r~*J 

_ , . JV 2 -1 

(2vrt)^— 

The factor of |j comes from the measure factor D/J, but for the limit of small t and finite N 
that we are considering for the moment we can drop this factor. Thus the result is that the 
small t limit of the diffusion probability is the kernel (see Appendix A) 

p(t,9)~lC(t,H) , (111) 

which is the fundamental solution of the heat equation on the full Lie algebra su(N). But as 
explained in Appendix A this is the free Brownian motion probability on a space of dimension 
N 2 - 1: 

N 2 -l 

K{t,e i )[dH] = k N 2_ 1 (t,9 i )de 1 ...de N 2_ 1 = ]J p^MOd* . (112) 

i=i 

Thus we recover the free Brownian motion on the tangent space of the group, which is what we 
expect. The second form in (112) is the product of the diffusion probabilities of N 2 — 1 abelian 
[7(1) theories. The Wilson loop distribution in that theory was computed in [13]: 

p u{1) (t,0) = -=e-£i (113) 

V 47TK 

where 

k = — =- , (114) 
4 v / 2^e 

and e is the smearing of the loop, which acts as a regulator. Thus, we get that 

t=^,0-0. (115) 
The effective action in this regime is thus 
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Let us comment on the limit t — > that we took above. It was useful to do the Poisson 
resummation (107) because it gives us an expression where we only need to keep a single term 
in this limit. If nevertheless we wanted to work directly with (106), we could have replaced the 
sum over representations by an integral. This is readily done by rescaling Aj by \fl in the usual 
way, so we consider large values of \ and keep A^ = \i\Jt fixed. Notice that in doing so we are 
not taking any large N limit. The integral one then gets is the usual Hermitian matrix model 
with Vandermonde interaction. Alternately, one can deform the theory by considering p(t; 9, 9') 
with non- vanishing 9'. One can then apply the techniques of [20] to get a Stieltjes-Wigert matrix 
model. At the end one can take the 9' — > limit. 

6.3 Casimir scaling regime 

Next we look at the Casimir scaling regime for large loops, where the potential is still irrelevant 
and we have an area law. The Casimir in the fundamental representation fj is C^fJ) = N—l/N, 
so we find that 



for large p. Thus, assuming that like in [13] that the logarithm of the Wilson loop is a sum of a 
constant, perimeter and area terms, the shape of the function r{p) that respects the two limits 



7 Vortex density and links with Lattice Gauge theory 

Once the spectral density is known one can also consider the simple vortex density; this is a 
quite interesting quantity to compute considering the importance of the center symmetry in the 
confinement problem and for our set up as well. In [11, 12] the SU (2) case for Wilson loops in the 
fundamental representation was addressed. For simplicity we consider here the same situation. 

Recall that fluctuations of the gauge group center, namely Z 2 in this case, are expected to 
be suppressed as soon as one moves to weak coupling but becomes pretty relevant at strong 
coupling. In particular various Monte Carlo simulations [12] show an almost exponential fall off 
(we are going to comment on this just below) for j3 > 2 of the thin vortex density E. The latter 
is defined as follows 



(115) and (117) is 




(118) 



E = \{1~ (sign(^))) 



(119) 



where j = 1/2 of course for the fundamental rep we are interested. In our case then Wj =1 / 2 
cos and one easily gets 




(120) 



Inserting our expression for G(9,t) one gets 




(121) 



n=2 
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The first term in our results is exactly the same found in [12]. We get qualitative agreement 
for the plot of E vs (3., i.e. the density decreases as soon as (3 grows as expected. It is quite 
difficult however to make a quantitative agreement for the following reasons which have to do 
with dimensional reduction. First of all, the four dimensional (3 and the two dimensional fod 
are not the same but schematically (3 = bfod as said. For small (3, b — > 1 but for large (3 then 
b ^ 2. This, as explained, is the way in which dimensional reduction works at the level of the 
couplings. Therefore in the Monte Carlo simulations one can only fit the data -after- choosing 
b in a suitable way. In our case we have an additional issue. Namely, the time t is inversely 
proportional to [3 but also depends on the size of the loop. This is additional residuum of the 
four dimensional physics in the dimensional reduction scenario. Therefore it sounds quite hard 
to provide quantitative matching with Monte Carlo data, but at qualitative level our model 
provides, consistently, the same leading order value at strong coupling and a progressive fall off 
as soon as one goes to weak coupling. 

Links with lattice gauge theory 

We believe that the lattice gauge theory and in particular the geometrical LGT approach can 
give at this stage a complementary and corroborative perspective to the whole scenario. This is 
also useful to explain a little better the meaning of the time t entering in the diffusion process. 

Geometrical LGT actions indeed ([21] and references therein) are supposed to describe fluctu- 
ations of plaquettes in the weak coupling regime; they are indeed weak coupling approximations 
of the standard LGT Wilson action. The prototype is the abelian Villan action [22] for magnets 
which is eventually generalized to a non abelian group as 



where the sum is over all irreducible unitary representations of dimension d r ; Xr is the character 
of the plaquette U, C r the (quadratic) Casimir in the representation r and N depends on the 
group we are considering. This action clearly satisfies the heat equation on the group manifold 
with initial condition S(U, 1) = exp(Sj(U, 1)). From physical point of view we imagine to study 
small (i.e. quadratic fluctuations) of the plaquette U around unity in group space; we use the 
geodesic distance of the plaquette from unity to measure them and the metric is induced on the 
group manifold by the invariant quadratic form of the Lie algebra. 

"Time" has therefore to be interpreted as the response to the heating and goes like 1/(3. 
Therefore the large time behavior we have been discussions means strong coupling regime and 
the initial conditions that one is starting from weak coupling to measure heat propagating. 

However, after we define the geometrical action we have to put it on a specific d-dimensional 
lattice configurations as usual in LGT to -concretely- compute. Therefore the total action will 
actually be the product of the single action one for each plaquette and the Haar measure entering 
in a path integral description will carry link indices. As usual it is clearly not possible to do an 
exact computation. If on the other hand one invokes the dimensional reduction hypothesis the 
geometrical action becomes precisely only at this stage the one of two dimensional QCD2 which 
is a character expansion indeed and one can forget now about all the plaquette links indices. 7 

We would like to make a simple example with the LGT geometrical action to reproduce 
the result of [13]. Consider free SU{2) for simplicity. The vev of the wilson loop in the spin j 

7 As an aside remark, note that the dimensional reduction is supposed to work not only for d — 4 — * d — 2 but 
for every d —* d — 2. [10]. 




(122) 
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representaion is 




(123) 



The geometrical action can be easily computed and becomes in the SU(2) case 



oo 



sin((2j + 1)0) 
sin(0) 



jjj + 1) 

2(3 



exp(- 



) 



(124) 



2j=0 



and 



1 



(125) 



(2j + 1) 



coincides with the SU(2) character in the spin j representation. Inserting into (123) one gets 



which is precisely what obtained in [13] at the qualitative level in the free diffusion case. 

8 Discussion and conclusions 

Building on the results of [13], we have presented a phenomenological model that describes 
the transition from Casimir scaling to screening. Motion on the group manifold SU(N) with an 
arbitrary potential breaks the simple Casimir scaling. The assumption of center symmetry of the 
potential predicts that the stable strings are the /c-strings. We have shown that asymptotically, 
the tension of the string does not depend on the representation of the Wilson loop, but only on 
its A-ality class. For intermediate areas, the expectation value of the Wilson is a sum over all 
possible representations of the same iV-ality, and this involves non-trivial dynamics. At large 
N, we expect the expectation value of the Wilson loop to be given by a two-matrix model. We 
will get back to this in the future [23]. It would also be interesting to see if phase transitions of 
the Douglas-Kazakov type are robust under perturbation by an arbitrary potential. 

A generalization to include time-dependent potentials as suggested in [13] should be straight- 
forward. One can then use time-dependent perturbation theory as in quantum mechanics. De- 
tails of this will appear elsewhere [23]. Our methods can also be applied when the potential no 
longer respects center symmetry; this is the case if dynamical quarks are introduced [13]. 

Our model is consistent with a closed string expansion in powers of N. Namely, it has an 
expansion in even powers of 1/N like the large N expansion of QCD2 [24, 25]. We want to 
stress however that our model is not QCD2. QCD2 is asymptotically trivial for large areas; the 
partition function goes to one, and there is no dynamics for large Wilson loops. Our theory 
however is non-trivial in this limit, see (85): the leading contribution to the Wilson loop is still 
interacting. 

One important question is what is the actual form of the potential? The simple symmetry 
requirement (48) was enough to obtain the right asymptotic behavior. To say more about the 
actual form of the potential, however, more phenomenological input is required. We leave this 
for the furure. Here we have shown that the model with general potential is solvable in group 
theoretic terms once the potential is provided. This is an interesting model to study in its own 
right. Providing the potential amounts to giving a finite set of numbers u a . The model in 



(Wj) = exp(- 



j{j + 1) 

2(3 



) 



(126) 
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principle contains enough freedom to accomodate for various values of the string tension. The 
only actual restriction is that the tension depends only on the A-ality class of the representation, 
and its leading contribution is given by the Casimir of the antisymmetric representation of that 
iV-ality. Because of this freedom, this model might in principle be able to mimic other kinds of 
behavior adequate to for example supersymmetric models [26, 27], where the string tension is 
given by the sine-formula: 



rather than the Casimir formula. As argued in [27], both formulas give the same result as 
N — > oo, namely k times the tension of the fundamental string (of A-ality one): 



Thus, the difference between both formulas is of order 1/N, and in principle this could be taken 
as the parameter that measures the perturbation of our potential. Thus, the phenomenological 
potential of [13] could in principle interpolate between both models. In fact, in one of the 
examples considered in Section (5.6) the correction to the ground-state energy came out to be 
the square of the quantum dimension, which applied to the tension of the stable string gives the 
quantum dimension of the antisymmetric string. This is already very reminiscent of the sine 
formula. However, we should remember that this is only a correction above the ground state 
energy. A more promising approach seems the following. Notice that both the sine law and 
Casimir scaling have invariance under k ^ N — k. If we view k as the level of some conformal 
field theory, this is level-rank duality. A second remark is that (127) is the quantum-deformed 
level, where the g-deformation parameter is q = e~ . Thus, it is tempting to speculate that 
the sine formula is somehow related to a ^-deformation. In fact, if we considered motion of the 
particle on the quantum group manifold, repeating the steps in this paper would replace the 
Casimir by the quantum Casimir (see for example [28]). 

There are many points that deserve further research. First of all, the expectation value of 
the Wilson loop computed in (83) itself, and not only the probability distribution, satisfies a 
heat equation. Presumably this heat equation can be seen as a loop equation. It would be 
interesting to analyze this in detail. 

It would be interesting to see if, apart from the phenomenological and mathematical interest 
of the model in [13], it can teach us something about dimensional reduction and the underlying 
theory describing screening. It is likely that the potential can be seen as an effective action 
that arises after integrating out some degrees of freedom. These degrees of freedom could be 
matter in an SU(N) representation of zero 7V-ality like the adjoint representation. Also, if 
we think of this matter as interacting with the gluons in a dimensionally reduced theory, in 
the screening phase the background necessarily will have to break the invariance under area- 
preserving diffeomorphisms. Simple cases of such models have been considered in [29]. It would 
be interesting to pursue this further. 

In addition it would be interesting to explore more in detail the Wilson loop effective action 
in the SU(N) along the lines sketched in Section 6 and in particular study in detail the large N 
limit of the model. 

Finally one could try to make a link along the lines of [30, 31, 32] with vicious random walk 
models trying to make a connection with the results of [33] . We leave all these topics for future 
research [23]. 




(127) 



a k = kai . 



(128) 
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A Laplacians and heat kernels on group manifolds 

In this Appendix we discuss the heat equation on the group manifold and the heat equation on 
the algebra(i.e. the full algebra and its Cartan subalgebra) we considered in the in the paper. 
The group manifold is in our case SU(N), namely the color group of the gauge theory. 

We first recall the expression for the radial Laplacian A on SU(N) which we considered in 
the paper 

A 1 1 9 J2 ( n ] 9 

2J 2 {9)d9i y 'd9i 
- ' ^-j { 9) + ±-N(N*-l) 



2J{9) dOf w 24 
J(9) = \[2sm±{9 i -9 j ) 

i<j 

m) = no*-*,-) 

i<j 

(129) 

where J (9) and D(9) are Jacobian the Vandermonde determinant of the Hermitian and unitary 
U (N) matrices. To get SU (N) one simplify the N variables 9i add up to zero: J2iLi @i = 0- 

Note the term proportional to N(N — 1) is just the curvature of the manifold and being a 
constant it can be shifted away. As noticed before we consider only class functions therefore 
we restrict to the radial Laplacian. For a rigorous derivation of the full Laplacian including 
the angular part see [34] and also [35] for a simplified treatment. We report in any case the 
expression of the angular part which is given by 



-E 

i<j 



4 + ^4 



or „ ( 13Q ) 



16sin 2 [(^-^)/2] 



with iLij = Eij — Eji and iMij = iiEij + Eji), where the Eij matrices are zero everywhere 
except one at the entries. The diagonal part is spanned by the iV — 1 Cartans while the 
Lij and My are sort of generalized iV 2 — iV step operators: in the simple SU(2) case they are 
just the usual J + and J_. 

We now recall some useful results from [21] concerning heat kernels. The fundamental 
solution of the heat equation 

(d t - A) K(t,9) = (131) 
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where A is the group Laplacian, is given by 



K(t,0) = £;aWXA(0)e-**W = 7 -^ I E D j^ + + |fj e-A(^) a . (132) 



A (27Tt) 

We also consider the solution 



(2vrf)^r- J ^ 



which is not, however, periodic [16] in 0. The relation between both is 

K(t,0i) = J2 K o(0i + ^h) ■ (134) 

k 

Thus, we see that Kq corresponds to the kernel in the covering space. In that sense, Kq lives in 
the algebra whereas K lives on the group. 
At small t, they are identical: 

hmK(t,<f>) = K Q {t,4>) . (135) 

The i-dependence above suggests that this solution is related to the solution of the heat 
equation on the space of Hermitian matrices: 

(d t -A H )lC(t,H)=0 (136) 

where 



-±TrH 2 



Klt,H) = - -TrT^e'-t 

V 1 (2^)^/2 

= I e -£££i«? (137) 

In the first form, H is a matrix and so the trace depends N 2 — 1 of its entries. In the second 
line we have done a similarity transformation = S*diag(^i, . . . , 9n)S* with a unitary matrix S. 
Notice that if U = exp(iH) is a unitary matrix, S is the same transformation that diagonalizes 
U: U = 5diag(e ie %...,e^)5t. 

It is now clear how both kernels are related: 

K (t,0) = jjQlC(t,e) (138) 



and 



At small t, 



Jr(»i) = 7^K(»i) = X"o(9i)- (MO) 



^Recall that in the SU(N) case J2i=o = °- 
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Finally, we define the Lie algebra kernel 

k N -!(t,0) = L-j-e-K^iO? (141) 

(2irt) — 

which satisfies the heat equation on the Cartan subalgebra of SU(N): 

($-A c )fc;v-i(t,0) = O. (142) 

The above is just the product of N one-dimensional Brownian motions. 
The kernels are correctly integrated as follows: 

1 = / [dH]lC(t,H) = [ dw N2 ~ 2 lC(t,<p) = [ d9 1 ...d9 N k N - 1 (t,9) 

JRN 2 - 1 JR Jrn- 1 

/N 2 -N f 
d0K(t,9) = (2irt)—*- [dU]K(t,U) (143) 

[dH] is as usual the Hermitian matrix measure, d9 the diagonal part of the unitary measure, 
d9 = d9\ . . . d9isfJ(6i), and [dU] the unitary measure [dU] = [dS] d9 where S is the similarity 
transformation used to diagonalize U. 9 denotes a radial coordinate for the N 2 — 1 variables of 
H. In particular, as densities we have 

K(t, <p) f N2 ^ d(p = k N -i(t, 9) d9 1 ...9 N , (144) 

which is what we used in the text. 

B Some elements of group theory 

We review some group theory notions which are used throughout the paper following the nota- 
tions and conventions of [36]. 

A weight written in the fundamental weight basis is: 

N N 

i=i i=i 

where are unit vectors in R N and 

K = iV E to ■ ( 146 ) 

j'=i 

So k = kN (mod N) is the iV-ality. The relation between both bases is: 

JV-l 
j=i 

i . N 

3=1 3=1 

for i = 1, . . . , N — 1, and ijy = 0. 1% is just the number of boxes in the ith row in the Young 
tableau of the representation A. Notice that, defining 

l i = e i -K + p i (148) 
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where the Weyl vector is 



N + 1 

ft = ~ * (149) 



we have 

N 



J2h = 0. (150) 



In the new basis, the inner product is now 

( A >^) =^2^iH( c ~ 1 )ij = ^2km . (151) 

ij i 

An SU(N) character is now: 

det ( e ilid - 



J ) - - 



»<*•■■ - , »>-dSF 5, 5r (152> 



and 

N N 



i=i i=i 

In this basis, the Casimir is given by: 

N 1 

C 2 (A)=^/, 2 --AT(Ar2-i). (154) 
i=i 

The Casimir in the trivial representation is zero. The Casimir in the fundamental is 

C 2 (D) = iV-l. (155) 

A commonly used convention differs from our by a factor of ^ . The Casimir for an antisymmetric 
representation of iV-ality k in our conventions is then 

C 2 (A(k)) = ^+1 k(N - k) . (156) 
In the main text we used the following representation of the delta function on a group: 

5(U) = ^djm(\)xx(U) (157) 

A 

where U is a group element SU (N) . The dimension is given by 

n ( *~(?- + o~° (i58 > 

l<i<j<N U ' 

and ii is the number of boxes in the ith row of the Young tableau. 
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It is easy to see that the center of SU (N) acts on the Cartan angles as $i like 6i — > 6i + ^ , 
for z = 1, . . . , TV. To see the effect on the character, we use the determinantal expression (152). 
We rewrite the determinant as: 

N 

det(0i) = det(e^)y = E e ( a ) II ( 159 ) 

treSjv i=l 

where 0jv = — (#i + ■ • • + #jv-i)- Shifting the determinant by the constant shifts gives: 
det(^ + — ) = Yl ' CT(1)+ ^ ' CT(2)+ - +2 ^ ^(^-D-^^-^W J] e a -(0ft 

creSjv i=l 
= e ( cr ) e 2 f i E£i' CTWe -2rf CT(iv) JJ e »I ff( i)9i ( 160 ) 

Notice that the first factor vanishes because X)£Li ^ = 0- The second factor gives: 



(161) 



for any i = 1, . . . , JV — 1. Therefore, we get 



det (^ + ^) = e 2 -^ e --( JV + 1 ) det(0i) (162) 

where = Nk is the iV-ality. We can now do the same in the denominator by setting k = 0. 
Thus: 

det(e^ + ^) = e-^W 1 ) det(e^) , (163) 
Therefore the minus signs cancel, and we are left with 

X x(0i + ^) = e 2 ^ N X x(9 i ) (164) 
which is what we wanted to prove. 

C Properties of the interacting theory 

We can use orthogonality of SU(N) characters of irreducible representations: 

(x», xu) = J ae j\e) xn(o)x,(o) = <W • (165) 

The inner product is defined with respect to the Haar measure, as usual. We now also require 
that the ip\S are orthonormal: 

<V>A,YV> = ^- (166) 
The same is true for the w's. This implies 

E 44= E <<> = ^ ( 167 ) 

M M 
k M- =k \ =k \> k M- = k \ = k \l 
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so that 

4 = ^ = (c-^, (168) 

so c is a unitary matrix. 

Having introduced g and n, a natural question to ask is: if we are free to choose the it's, how 
many components of n can we choose? Notice that 

< = u„ (169) 

where is the trivial representation. Therefore, choosing is equivalent to choosing Uq. Once 
this is done, all other components of riu are fixed. 

We can use the charge conjugation matrix to lower indices: 

c m p = . (170) 

Indeed, this definition is correct because 

JV"i = %ux = J [d9] x»(0)Xu(0)xx(e) ■ (171) 

It therefore satisfies 

Jvi = Kx • ( 172 ) 

From this, we derive 

»v = n ^ = n v ( 173 ) 

and 

< = 4 • (174) 
Notice that the newly defined eigenfunctions ijj also satisfy: 

A 

The TV's can be computed in terms of the tensor multiplicity coefficients: 

iV^ = 5>^XV • (176) 

ctct' 

Finally, we should ensure positive definiteness of the potential. This means 

{^\U\^) = J2^i'^>0. (177) 

AA'cr 

This ensures that the energy is positive. Indeed, from the above we get for the additional term 
in the energy: 

5E X > . (178) 
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